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METHOD OF SEGMENTING A RADIOGRAPHIC IMAGE INTO 
DIAGNOSTICALLY RELEVANT AND DIAGNOSTICALLY 
IRRELEVANT REGIONS 



5 FIELD OF THE INVENTION 

This invention relates in general to the processing of digital 
radiography, and in particular to the segmentation of a digital radiographic image 
into diagnostically relevant and diagnostically irrelevant regions. 

1 0 BACKGROUND OF THE INVENTION 

The recent advance in digital radiographic imaging systems, such 
as the flat-panel based direct digital radiographic (DR) systems ("A high 
resolution, high frame-rate, flat-panel TFT array for digital x-ray imaging/' 
Proceedings of SPIE Medical Imaging, Antonuk et al., vol. 2163, pp 1 18-128, 

15 1994) and the storage phosphor based computed radiographic (CR) systems 
("Introduction to medical radiographic imaging," Eastman Kodak Company, 
1993) etc, allow the separation and optimization of image acquisition, processing, 
and display processes individually. The image pixel data are manipulated during 
the image-processing step so as to optimize the image information acquired on the 

20 radiograph and to help the radiologists to better perceive even the subtlest 

diagnostic details. Optimizing diagnostic details depends on the knowledge of the 
location and characteristics of both diagnostically relevant and diagnostically 
irrelevant regions in the radiograph. The scope of this invention therefore relates 
to the automatic segmentation of a digital radiograph into anatomy (diagnostically 

25 relevant regions), foreground and background (diagnostic irrelevant regions). 

Fig. 1(a) shows an example of a foot radiograph acquired with CR. 
The foot was exposed at three viewing angles A, B, C using the same storage 
phosphor cassette but at different cassette regions. The anatomy in this 
radiograph is the foot and is considered diagnostically relevant. The regions in the 

30 radiograph where x-rays directly expose the storage phosphor are diagnostically 
irrelevant, which are later referred to as the background or direct exposure region. 
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Collimation was used during the x-ray exposure to reduce unnecessary radiation 
to the anatomy that is irrelevant to diagnosis and to confine the x-rays to a local 
region of the cassette. The regions in the radiograph collimated outside the x-ray 
radiation field are diagnostically irrelevant, which are later referred to as the 
5 foreground region. 

The work of Barski et al. taught a method of determining the direct 
exposure regions in a radiograph ("Determination of direct x-ray exposure regions 
in digital medical imaging," U.S. Patent 5,606,587). However, this method uses a 
single value, namely background left point, to threshold the input image to 

1 0 generate a binary map. Although this method is suitable for the purposes for 
which it was intended, in the case when the direct exposure region has a large 
intensity non-uniformity, which can be caused by the non-uniform intensity 
distribution of the incident x-rays or by the exam requirements that need multiple 
views of the anatomy at different exposures levels another method may be needed. 

15 Fig. 1(b) shows an example of the result using Barski 's method, which indicates 
that the background can be over-detected and under-detected at the same time. 
The thresholds used in this method also vary from one exam-type (bodypart and 
projection) to the other, which may not work well when exam-type information is 
unavailable. 

20 Pieter disclosed a method to determine the foreground region in a 

radiograph ("Method of recognizing an irradiation field," European Patent 
0,610,605). However, this method can only deal with single-exposed radiographs, 
and can not be used to completely segment the diagnostically relevant regions 
where there are two or more radiation fields in a single radiograph, such as the 

25 example in figure 1(a). To address this limitation, Piet et al disclosed an improved 
method of detecting the foreground regions ("Method of recognizing one or more 
irradiation fields," European Patent 0,742,536). However, both Pieter and Piet 
failed to teach a method for determining the background regions. Wang et al., 
showed a method based on Barski 's direct exposure detection to determine the 

30 image foreground regions ("Method for recognizing multiple irradiation fields in 
digital radiography," US Patent 6,212,291). This method, however, still did not 
fulfill all of the needs of background detection. 
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Given the drawbacks and limitation of the prior art, there is a need 
for a method that is exam type independent and that can automatically segment 
digital radiographic images into anatomy, foreground, and background regions. 

5 SUMMARY OF THE INVENTION 

According to the present invention, there is provided a solution to 
the problems and a fulfillment of these needs. 

According to a feature of the present invention, there is provided a 
method of segmenting a radiographic image into diagnostically relevant and 
10 diagnostically irrelevant regions comprising: 

acquiring a digital radiographic image including a matrix of rows 
and columns of pixels; 

detecting the initial background left point of a histogram of said 

image; 

1 5 detecting the foreground of said image; 

regenerating the background of said image by region growing; 
validating the background of said image; 

merging the background and foreground regions of said image as 
diagnostically irrelevant regions; and 
20 extracting the anatomy region of said image as the diagnostically 

relevant region. 

ADVANTAGEOUS EFFECT OF THE INVENTION 

The invention has the following advantages. 
25 1 . It is independent of the exam type of the radiograph. 

2. It can automatically segment a radiograph into foreground, 
background and anatomical regions with no over- or under- detection. 

3. It can successfully detects the whole image background 
using region growing regardless of the background non-uniformity. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 a is a diagrammatic view which shows an example of a foot 
radiograph acquired with CR, and the definitions of image foreground, 
background (direct exposure), and anatomy. 
5 Fig. lb is a diagrammatic view which indicates over- and under- 

detection of image background using a known method. 

Fig. 2 is a flow chart of the major processing steps in one 
embodiment of the present invention. 

Fig. 3 is a flow chart of the major processing steps in an 
1 0 embodiment of the present invention for background detection. 

Fig. 4a is a diagrammatic view which uses a dotted line to 
represent a sample of an image row. 

Fig. 4b is a graphical view of the profile of the sampled line of Fig. 
4a (darker image regions have higher pixel values). A transition sample is marked 
15 by the ellipse. 

Fig. 4c is a graphical view of a modified presentation of the ellipse 
region of Fig. 4b. Some definitions, such as maxTransMag, segLen, segRange, 
high _j?t, and maxPixelValue y are also indicated. 

Fig. 5 a is a diagrammatic view which shows an identified 
20 candidate background pixels. 

Fig. 5b is a graphical view which shows an image global histogram 
and a histogram of the candidate background pixels. 

Fig. 5 c is a graphical view which magnifies the candidate 
background pixel histogram of Fig. 5b and defines peak,peak_lp,peak_rp, and 
25 last jft. 

Fig. 6a is a diagrammatic view of an example of the gradient 
magnitude image. 

Fig. 6b is a graphical view of the histogram of the Fig. 6a image. 
Two parameters, peak_hw, and peak_loc, are defined. 
30 Fig. 7 is a graphical view showing an image global histogram, in 

which pValuemin, pValue m ax, and bklp__lp are defined. 
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Fig. 8a is a diagrammatic view which shows the detected initial 
background regions. 

Fig. 8b is a diagrammatic view which show foreground regions 
detected based on Fig. 8a. 
5 Fig. 9a is an image with foreground removed (the darkest region) 

then overlaid with the initial seeds for background region growing (the dark dots 
in the background regions). 

Figs. 9b and 9c are diagrammatic views which show some 
intermediate results of background region growing. 
10 Fig. 9d is a diagrammatic view which shows the final detected new 

background regions. 

Figs. 10a and 10b are diagrammatic views which show examples of 
the detected background and foreground regions, respectively. 

Fig. 10c is a diagrammatic view which shows the sum of 
15 foreground and background, from which the transition gaps are seen. 

Fig. lOd is a diagrammatic view which shows the image with 
background and foreground merged together. 

Fig. 1 1 is a diagrammatic view which shows the final anatomic 
region detected after the merged background and foreground are removed. 

20 

DETAILED DESCRIPTION OF THE INVENTION 

Referring now Fig. 2, an embodiment of the present invention will 

be described. 

A digital radiographic (medical) image is acquired (box 201) such 
25 as from a diagnostic imaging device (MRI, CT, PET, US), from a direct digital or 
computed radiography device, from an x-ray film digitizer, from a medical image 
archive. The digital medical image includes an array or matrix of rows (lines) and 
columns of pixels having a gray scale range of a predetermined range of bits or 
code values (e.g., 8 or 12 bits). The digital medical image is processed in an 
30 image processor, such as a digital computer and output to a display or hard copy. 
The method of the invention is carried out in the image processor. 
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This invention first tries to improve the background detection 
algorithm disclosed in U.S. Patent 5,606,587, issued February 25, 1997, inventors 
Barksi et al., by making it exam-type independent and at the same time to provide 
a reasonable initial background left point (box 202) for the acquired digital 
5 radiographic image. Then it takes this initial background left point to threshold 
the digital radiographic image in order to create an initial background map, which 
is used as one of the input for applying the foreground detection algorithm 
disclosed in U.S. Patent 6, 212,291, issued April 3, 2001, inventors Wang et al. 
(box 203). This results in a well-detected foreground region and a less-than- 

10 desirable background region. Later on, a region growing process is performed 
from some identified "seeds" to regenerate the whole background map (box 204). 
After the new background is validated (box 205), the background and foreground 
regions are merged together (box 206) and finally the rest of image region is 
identified as anatomy (box 207). 

1 5 The flow chart of the initial background detection algorithm is 

shown in Fig. 3. In summary, this part of the algorithm tries to find a good sample 
of image background pixels, from which a reasonable value (later referred to as 
initial background left-point) is determined which can be used to binarize the 
image into background and non-background regions. The background determined 

20 may not be perfect, but is sufficient for the foreground detection step (box 203). 

The process starts with analyzing each column and row of image 
pixels for all the rising and falling transitions (box 301) in order to find all 
candidate background to foreground (BF) and background to skinline (BS) 
transitions. The horizontal dotted line 400 in Fig. 4a represents a sample row of 

25 image pixels, where those pixels in the darker region have higher pixel values. 
The profile of the sampled row is plotted as line 402 in Fig. 4b, in which a rising 
transition sample is marked by the dotted ellipse 404 and magnified in Fig. 4c. As 
shown in Fig. 4c, there are four major measurements associated to each transition, 
transition width (segLeri), transition range (segRange), maximum transition slope 

30 (maxTransMag), and transition high-point value (high _pi). Each of the four 

measurements is computed and recorded individually for all the transitions in the 
image. The BF and BS transitions tend to have a larger transition length, 
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transition slope, and maximum transition magnitude, plus the transition high-point 
value tends to have a larger value. Therefore, the following criteria are applied to 
eliminate those false candidate transitions that may correspond to bone-tissue, 
tissue-tissue, and any other transitions: 
5 segLen > threshold__segLen, 

segRange > threshold_segRange y 
maxTransMag > threshold_JransMag y 
high _j)t > maxPixelValue - A> 
where maxPixelValue is the maximum pixel value of the image, and A is a 

10 predetermined value which ensures that it is large enough to include all the BF 
and BS transitions, and at the same time small enough to exclude as many as 
possible any other transitions. 

The aforementioned thresholds are generically created for all 
exam-types. They can be optimized accordingly when exam type information is 

1 5 available. However, this invention does not have to rely on exam-type 

information, it extracts a set of features from the image itself then adaptively 
detects the image background. 

The resultant transition candidates in this processing step may still 
contain some false detections other than the valid BF and BS transitions. Because 

20 the segRange values of the undesired bone-tissue and tissue-tissue, etc. transitions 
are relatively small compared to those of the valid transitions, they will most 
likely be distributed at the low end of the segRange population. Based on this, the 
histogram of segRange is built, from which the cumulative histogram is 
computed, and a new threshold segRange is set at a certain percentile such as 50% 

25 of the cumulative histogram (box 302). This new threshold is used to further 

prune the transition candidates. Slightly over-pruning is not a problem as long as 
there are enough candidate transitions left for subsequent processing. Because 
images of different patients and different exam types usually have different BS 
transition characteristics, the selection of a new threshold segRange based on a 

30 certain percentile of the cumulative histogram can automatically adapt the pruning 
process based on the image characteristics itself, i.e., making the algorithm 
independent of exam type information. 



-8- 

The high _pt of the detected candidate BF and BS transitions are 
considered to be candidate background pixels in the image. However, the 
background region for some images may encompass an entire image line (row or 
column) and there is no significant transition that can be found within the line 
5 using the aforementioned transition detection process. To solve this problem, 
first, a pixel value threshold (thresh_column and thresh_j-ow) is defined for each 
image row and column: 

thresh_column [i] = minimum high _pt of all candidate BF and BS 
transitions in column /, 
10 thresh_row [/] = minimum high _pt of all candidate BF and BS transitions 

in row j. 

Second, those image pixels whose values are greater than both their 
corresponding row and column thresholds are considered candidate background 
pixels. An example of the candidate background pixels is shown in Fig. 5a as 

15 coarsely distributed white dots 500. The image is considered to have no 

background if there are not enough number of candidate background pixels are 
found at this processing step. 

A histogram of the high _pt for the identified candidate background 
pixels is built for the estimation of the overall characteristics of the image 

20 background (box 303). Because the background pixels have relatively higher 
values among all the image pixels, in a preferred embodiment of this invention, 
only those candidate background pixels whose high jpt values are within a certain 
range below the maximum pixel value of the whole image are included in the 
histogram calculation. Fig. 5b shows an example of the candidate background 

25 pixel histogram overlaid on top of the image global histogram. Fig. 5c shows the 
magnified candidate background pixel histogram, from which four features are 
extracted: peak, peakjp, peak_rp, and last _j?t. Feature peak is the location of the 
histogram peak, peakJLp (peak_rp) is the peak left (right) point, which is defined 
as the location in the histogram whose population exceeds a certain threshold, and 

30 last _j?t is the last (highest) pixel value in the histogram. It needs to mention that 
these four aforementioned features will be used in the iteration process ( box 306). 
If last_pt is too far below the maximum pixel value of the whole image or if the 
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distance between last _jpt and peak_jrp is too far, it is considered that background 
does not exit. Otherwise, the resultant peakjp value is taken as the first estimate 
of the background left point and is used as the threshold to get a first estimate of 
the image background region. Any image pixel having a higher value than 
5 peak_lp is considered as a background pixel, and all the rest of the pixels are 

considered as non-background pixels. The true background region, by definition, 
should be relatively smooth because this is the region where the direct x-rays 
expose the image receptor. The pixel values in the true background region 
therefore should be at the higher end of the image histogram and their distribution 

10 tends to be within a small range. Observation of the image histogram supports 
this judgement (Fig. 5b): the background pixels tend to create a sharp peak in the 
very high end of the histogram. This sharp peak is detected in the algorithm (box 
305) as one of the unique features: major _j?eak, where major j>eak = true if at 
least one sharp peak can be detected from the high end of the image histogram, 

15 and major jpeak = false otherwise. Feature major _j?eak is used in the iteration 
process (box 306) to select different thresholds. 

The distribution of the background pixel values is within a 
relatively small range suggests that their variation should be relatively small too. 
Therefore, if the gradient is calculated from the original image (box 304) then the 

20 background region should have relative small gradient magnitude (Fig. 6a). When 
the histogram of the gradient magnitude is plotted (Fig. 6b), a relatively narrow 
ipeakjiw) peak can be found near (peak_Joc) the histogram origin. Feature 
peak_loc itself represents the most populated gradient magnitude, and peak_Jtw 
represent how small the variation of the most populated gradient magnitude is. 

25 Based on peak_loc and peak_Jiw a new feature is defined (box 305): major_J>kgr, 
where majorjbkgr = true if peak_loc < thresh __peak_Joc and peak_Jiw < 
thresh_peak_hw, and major_bkgr = false otherwise. When majorjbkgr = true the 
image most likely has a significant number of background pixels present, 
therefore, this feature is also used in the iteration process (box 306) to select 

30 different thresholds. There are many filters that can be used to calculate the image 
gradient, such as Robert, Sobel, Previtt, or isotropic filters. In a preferred 
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embodiment of this invention, the gradient is calculated using the isotropic filter 
and square root is taken for the gradient magnitude. 

The candidate background left point is initially set being equal to 
peak__lp then fine-tuned in the iteration step (box 306). The iteration concept is 
5 one of the most important features of the present invention in making the 

segmentation independent of exam type information. In this processing step, any 
image region having a higher pixel value than the background left point is 
considered as candidate image background and is validated based on the regional 
gradient magnitude. In particular, a measurement is defined: 

1 0 bkgr^yariance = ^ gMag(i, j ) / ^ , 

tj u 

where gMag(i, j) is the gradient magnitude at pixel location (/, j) and the sum is 
over only the candidate background regions. Because the BF and BS transitions 
have fairly large gradient magnitudes, if the candidate background left point is 
defined too aggressively, i.e., not only the true background regions but also some 

1 5 regions of foreground or anatomy are included, the large gradient from BF and BS 
transitions can make bkgr_yariance really large. 

Based on this, the iteration starts with the initial background left 
point value ipeak_Jp)^ calculates bkgr_yariance and checks if it is greater than a 
certain threshold. If it does, the background left point value will be increased by a 

20 certain amount, and then the iteration restarts again. Every time the background 
left point value increases, the number of candidate background pixels decrease by 
a certain amount and so does the size of candidate background regions. The 
criteria used for stopping the iteration are: number of iterations, number of 
candidate background pixels, and value of bkgr__variance. These criteria depend 

25 on the aforementioned two features: major _j>eak and majprjbkgr so that the 

iteration can be controlled differently when either of major _peak or majprjbkgr is 
true. The iteration stops if there are not enough candidate background pixels left 
or if there are too many iterations. This means either that the image background 
most likely does not exist, or if the criteria for bkgr^yariance is met, that a 

30 reasonable background left point has been found. The decrement amount for 

background left point can be fixed in pixel value, can be a calculated value fixed 
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in the amount of candidate background pixel decrement, or can be a combination 
of both. In a preferred embodiment of the present invention, the combination 
method is used to ensure that the decrement is fine enough for finding the most 
likely background left point but still large enough to minimize the number of 
5 iterations. 

The final processing step is to verify the initial background left 
point (box 307). Fig. 7 shows an example of the whole image histogram, in which 
pValue min and pValue max are the minimum and maximum pixel values of the 
original image, and bkgr_lp is the background left point obtained from the 

10 iteration step. In a preferred embodiment of this invention, pValue min and 
pValue max are calculated from the cumulative histogram of the original image 
pixel values at some percentile levels, such as 0.5% and 99.9% for pValue min and 
pValue maX y respectively. One implementation of the verification step uses a rule- 
based reasoning process. For example, a valid background left point must satisfy 

1 5 the following criteria: (1 ) the value of bkgrjp should be sufficiently large, (2) the 
dynamics range of the image, as defined by pValue max - pValue min should be 
sufficiently large, (3) the area of the background regions should be sufficiently 
large if either majorjbkgr or major jeak is true etc, or should be sufficiently 
small otherwise. 

20 The first rule is used because radiographs are exposed at a certain 

predetermined condition and the direct exposure regions tend to have a relatively 
high exposure. The dynamic range of a radiograph depends on the exposure 
technique, such as kVp etc, and the human anatomy imaged. Because the 
anatomy introduces quite a bit of x-ray attenuation, the acquired radiograph 

25 should have a reasonably large dynamic range if there is any background regions 
in the image, which is the rational for the second rule. The reason for the third 
rule is because the image tends to have a large background region when either 
majorjbkgr or major__peak is true, or tends to have a small background region or 
none at all, otherwise. Finally, the number of candidate BF and BS transitions, 

30 whose high_pt values are higher than the initial background left point, is 

compared to the total number of candidate BF and BS candidates, and, if the ratio 
is too small, then the identified initial background left point is considered invalid. 
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Fig. 8a shows an example of the detected background pixels, which is not perfect 
but is good enough for the subsequent foreground detection process. 

The foreground detection processing used step is disclosed in 
"Method for recognizing multiple irradiation fields in digital radiography," U.S. 
5 Patent 6,2 1 2,29 1 , inventors Wang et al. In summary, the algorithm ( 1 ) takes the 
original image and the initial background left point value as inputs, (2) uses a 
smart edge detection process to classify all significant transitions in the image, (3) 
conducts Hough Transform to delineate all the lines that are possible collimation 
blades, (4) finds candidate partition blade pairs if the image has several radiation 

10 fields, (5) uses a divide-and-conquer process to partition the image into sub- 
images containing only one radiation field, (6) identifies the best collimation for 
each sub-image. Fig. 8b shows an example of the detected foreground. 

The background detected so far has many drawbacks, such as over- 
and under- detection as shown in figure 1(b) when the background region has a 

15 large non-uniformity. This processing step tries to overcome these drawbacks by 
taking advantage of the known foreground and regenerating the background 
region using region growing. In one embodiment of this invention, the new 
background region is regenerated using the region growing method from a set of 
known "seeds." These seeds can be the initial background pixels that have been 

20 detected before, or can be the detected background transition pixels etc. In one 
embodiment of the present invention, the criteria for region growing are based on 
the homogeneity between the seeds and the detected background region. It 
assumes that the pixel value distribution in the background region can be 
approximated by a Gaussian distribution. The method first computes an average 

25 (jiBackgound) and a standard deviation (&Backgowud of the detected background, then 
starting the region growing using a 3x3 neighborhood region of a seed. The merge 
condition is defined as follows: 

\M M Background j — 
^ 0 ^Background 

where the constant t represents a confidence interval in a normal distribution. If 
an average pixel value (ju) of a seed region falls within a predetermined 
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confidence interval, say 90%, the seed region will be merged into the background, 
otherwise, it is marked as anatomy or foreground region. Once all seeds have 
been checked, a set of new seeds are generated around the boundary of the newly 
segmented background. Growing of the seed regions continues until no region 
5 can be merged with the background. 

In another embodiment of the present invention, the detected 
candidate background transition pixels whose highjpt values are greater than the 
initial background left point are taken as the seeds. The region growing is based 
on the local pixel value gradient. In particular, the gradient is calculated based on 
10 a 2D detrending process. The original image is first partitioned into many 

partially overlapping 3x3 blocks, then each block is fit to a 2D bilinear function: 

bfayj) =axi + byj+ c, 
where xi and y } are the pixel coordinates, ie [1,3], je [1,3], and a, b, and c are the 
fitting parameters. Least square error 

15 A 2 = £[/(x J; >6(x J; .)] 2 

is used to find the best set of fitting parameters, where i(x if yj) is the image pixel 
value at (x it yj). The gradient magnitude for the center pixel of each block is 
defined as: 

grad = yja 2 + b 2 . 

20 The growing processing first starts with the known seeds, it then 

grows to its neighborhood if grad is smaller than a certain threshold. The rational 
is that the grad value should be relatively small in the background regions, but 
should be very large in the region near BF and BS transitions. The growing 
process merges all the background regions but will stop near BF and BS 

25 transitions therefore will not get into the anatomy region. 

An example of the region growing results are shown in Fig. 9a to 
9d. Fig. 9a shows the seeds 90 overlays on top of the original image with 
foreground removed. Fig. 9b to 9c show some intermediate steps as to how the 
candidate background regions 92 grows. The final result in Fig. 9d shows no 

30 background over- and under detection (as in Fig. lb). 
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There is a final validation step to check the background region 
grown from the seeds. Because the foreground is known and therefore can be 
eliminated from the analysis, a more accurate testing of the background pixels can 
be made. In particular, because the anatomy should be of a minimum size when 
5 background exits, the ratio between the background area and all the non- 
foreground area should not exceed a certain threshold, and the dynamic range of 
all the non-foreground area should be greater than a minimum threshold. 
Otherwise, the background detected is considered as invalid. 

Both the background and foreground regions in the image are 

1 0 diagnostically irrelevant. The diagnostically relevant anatomical regions are 

extracted by merging and then excluding the foreground and background regions 
from the original image. Figs. 10a and 10b show an example of the detected 
background and foreground regions, respectively, and figure 10c is the sum of the 
two. One can find that the foreground and background regions may not connect to 

15 each other so that there are transition gaps between them. These transition gaps 
must be removed in order to obtain a "clean" anatomical region. To do this, each 
column and row of pixels of the original image is analyzed. If a transition is 
found that begins with foreground then ends with background, or begins with 
background then ends with foreground, it is recorded. If within this transition 

20 there is a small number of pixels that are neither inside the background nor inside 
the foreground and are in the proximity of the transition maximum, this transition 
is used to "bridge" the foreground and background region. This process is 
repeated for every transition in each image row and column until all the transition 
gaps are removed. Fig. 1 Od shows an example of the resultant image. 

25 The anatomy area is extracted by subtracting the merged 

foreground and background regions from the original image. There may be some 
residuals left in the detected anatomy image that are diagnostically irrelevant, such 
as hardware markers etc. They are usually very small in size and may not 
constitute a problem for subsequent image processing. Fig. 1 1 shows an example 

30 of the detected anatomy region. In case hardware markers constitute a problem, 
the anatomy image can be further processed so as to refine the anatomical regions. 
For example, the anatomy image can be labeled such that each connected region is 
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assigned a unique pixel value. Because the anatomy should generally connected 
and should be larger than any hardware markers in size, a histogram of the labeled 
image is calculated, the peak in the histogram will identify the pixel value of the 
largest connected region in the labeled image, therefore the largest connected 
5 region can be extracted while leaving out those small and diagnostically 

unimportant connected regions (e.g., hardware markers). If there are two or more 
connected regions of similar size, one can extract all as diagnostically important 
regions. 

The invention has been described in detail with particular reference 
10 to certain preferred embodiments thereof, but it will be understood that variations 
and modifications can be effected within the spirit and scope of the invention. 
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90 seeds 

92 background regions 

201 processing step for acquiring digital radiographic image. 

202 processing step for detecting initial background left point. 

203 processing step for detecting image foreground. 

204 processing step for regenerating background by region growing. 

205 processing step for validating background. 

206 processing step for merging background and foreground. 

207 processing step for extracting anatomy region. 

301 processing step for evaluating significant transitions and recording 
segRange and high jpt. 

302 processing step for building segRange cumulative histogram and auto- 
selecting new segRange threshold. 

303 processing step for building high _pt histogram and finding first estimate of 
background left point. 

304 processing step for building image gradient map. 

305 processing step for identifying major _pedk & majorjbkgr. 

306 processing step for iteratively fine tuning background left point. 

307 processing step for verifying background left point. 
402 sampled row is plotted 

404 dotted ellipse 

500 distributed white dots 



